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We perform electrical field effect measurements at 4 K on insulating granular aluminium thin 
films. When the samples size is reduced below ~ lOO^tm, reproducible and stable conductance 
fluctuations are seen as a function of the gate voltage. Our results suggest that these fluctuations 
reflect the incomplete self-averaging of largely distributed microscopic resistances. We also study 
the anomalous field effect (conductance dip) already known to exist in large samples and its slow 
conductance relaxation in the presence of the conductance fluctuations. Within our measurements 
accuracy, the two phenomena appear to be independent of each other, like two additive contributions 
to the conductance. We discuss the possible physical meaning of this independence and in particular 
whether or not this observation is in favor of an electron glass interpretation of slow conductance 
anomaly relaxations. 

PACS numbers: 72.80.Ng, 61.20.Lc, 73.23.-b, 73.23.Hk 



I. INTRODUCTION 

An anomalous field effect in disordered insulators has 
been first observed in 1984 by Adkins et al. pQ. In 
their experiment, an insulating discontinuous gold film 
was used as the conducting channel of a MOSFET struc- 
ture. At low temperature, the conductance of the film 
was found to increase when the gate voltage was swept 
away from its cooling value, whatever the sweep direc- 
tion. A few years later, Ovadyahu et al. observed a sim- 
ilar conductance dip in insulating indium-oxide (amor- 
phous InOa; and crystalline I^O^x) thin films [2]. A 
large number of experiments have been performed since 
then in the latter systems [3]. The conductance was 
shown to decrease as a logarithm of time after a tem- 
perature quench or a sudden gate voltage change. The 
samples never reach a true equilibrium state even after 
one week of measurement. Qualitatively similar effects 
were reported in ultra thin films of Pb and Bi [3] and we 
published a thorough study of these phenomena in gran- 
ular Al films [5j[6]. It is worth noting that doped semi- 
conductors do not show a dip and that slow and glassy 
like electronic relaxations have been only rarely reported 
in there [3 [8]. 

A priori all these systems have different microstruc- 
tures: the disorder is "granular" for discontinuous 
and granular metal films whereas it is "homogeneous" 
(oxygen vacancies and atomic position distribution) for 
indium-oxide films. Consequently, two different physical 
pictures have naturally emerged to explain the anoma- 
lous field effect and its long relaxation times. 

- In InOa; films, it was shown jS] that the typical 
relaxation time and the width of the conductance 
dip depend systematically on the carrier concentra- 
tion n: the larger n, the longer the relaxation time 
and the wider the dip. Moreover, the relaxation 
time decreases rather sharply below n ~ 10 20 cm~ 3 . 
These results may explain why the anomalous field 



effect and its slow relaxation are observed in An- 
derson insulators with a high carrier concentration, 
like indium-oxide and granular Al, and not in doped 
semiconductors which usually have much smaller n 
[9] . They also suggest that such phenomena reflect 
the existence of an electron glass at low temper- 
ature (we call this hypothesis the intrinsic one). 
The electron glass existence was theoretically pre- 
dicted for disordered insulators in 1982 by different 
authors [THl EJ HI] • The combined effects of disor- 
der and unscreened interactions between localized 
electrons should give rise to a glassy dynamics and 
correlated motion of electrons at low temperature. 
The experimental findings of slow conductance re- 
laxations have motivated several theoretical works 
on the electron glass problem [T3 | [T4l H5 l [TB I [17] [18] . 

- In granular systems, a different explanation was 
proposed. It was shown that the conductance dip 
and its slow relaxation could be qualitatively ac- 
counted for by a slow polarization of the dielectric 
material around the metallic grains 6 . This idea 
was first introduced by Adkins et al. [T] and de- 
veloped later by Cavicchi and Silsbee jT^] for the 
interpretation of capacitive measurements on gran- 
ular films. We call this hypothesis the extrinsic 
one since according to it, the non equilibrium ef- 
fects don't come from the electrons of the metallic 
grains themselves. Other observations may be in fa- 
vor of an extrinsic interpretation. Features qualita- 
tively similar to field-gated conductance measure- 
ments were found at very low temperature in the 
dielectric response of amorphous materials to large 
DC electric fields [20]. It was also shown that bias 
voltage changes can induce at room temperature 
non exponential relaxations and memory effects in 
the electronic properties of Al-AlO^-Al planar tun- 
nel junctions, the effects being attributed to the 
metastability of interface states [2T]. Interestingly 
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enough, such tunnel junctions are present in our 
granular Al thin films [6|. Last, we would like to 
mention that an extrinsic type of scenario has been 
recently suggested for indium-oxide films |22J. 

The striking similarity between the different experimen- 
tal results strongly suggests that the same physical pro- 
cess is involved for all the systems. We have already 
discussed in Ref. [B] why we believe that the two hy- 
potheses (extrinsic versus intrinsic) could be applied to 
both granular and homogeneous systems. Up to now and 
to our point of view, none of the hypothesis has been un- 
doubtedly settled by the experiments. The logarithmic 
time dependence of the conductance and the simple aging 
relaxation laws observed after gate voltage changes can 
be interpreted as the response of independent degrees of 
freedom with relaxation times Tj such as In Tj has a flat 
distribution over the time scale of the experiments [6 , 23J. 
Such an interpretation is obviously independent of the 
nature of these degrees of freedom (correlated electrons 
hops, two level systems charge configurations, etc. |H]). 
One important difficulty in order to rule out or to estab- 
lish one hypothesis comes from the absence of quantita- 
tive and indisputable predictions directly comparable to 
experimental data. Whether the intrinsic (electron glass) 
hypothesis is the correct one or not is of prime interest 
since the long response times measured (more than hours 
and days) may constitute a direct experimental evidence 
of its existence. 

To shed light on the actual origin of these slow relax- 
ation phenomena, we choose to perform field effects mea- 
surements on small size granular Al thin films. By small 
size, we mean small enough to observe reproducible gate 
voltage induced conductance fluctuations and to study if 
the conductance dip is modified by these fluctuations. In 
diffusive metals of mesoscopic size, the conductance was 
shown to be sensitive to details of disorder and even to 
the motion of single impurities (see Ref. [25] for a review) . 
Changing the gate voltage or the magnetic field modifies 
the quantum interference effects between elastic diffusion 
centers responsible for weak localization corrections. At 
low temperature, this gives rise to the well known Uni- 
versal Conductance Fluctuations, which are of the order 
of the quantum of conductance e 2 /h for samples having 
sizes equal or smaller than the coherence length. Gate 
voltage or magnetic field conductance fluctuations pat- 
tern could therefore be considered as a finger-print of the 
sample specific realization of disorder. 

A similar finger-print exists also at low temperature for 
hopping disordered systems. Due to the large (exponen- 
tial) distribution of electron hops probabilities between 
pairs of localized states, the conductance of small size 
samples is dominated by a few critical hops or micro- 
scopic resistances The gate voltage and the mag- 
netic field could change the values and the positions of 
the dominant resistances, giving rise to conductance fluc- 
tuations of quantum |27[ 128] and/or geometrical origin 
[2"9"] . The typical distance between the dominant resis- 
tances defines the length scale above which the conduc- 



tance starts to self-average [26, 28 . This length scale 
could be orders of magnitude larger than the phase co- 
herence length (the hopping length in the variable range 
hopping regime) when the spread of hops resistances is 
very large. Indeed, conductance fluctuations have been 
observed in samples of sub-micrometer [30j |3TJ [32] mi- 
crometer [33J [31] and millimeter sizes [351 Tele- 
graphic noise of individual fluctuators was also identified 
[32l 1331 137] , Depending on the system parameters and 
the temperature, the amplitude of the fluctuations can 
be as large as the conductance itself. As we will discuss 
in more details in Sect. |III B] similar gate voltage conduc- 
tance fluctuations are expected in our insulating granular 
Al films. Since the fluctuations pattern obtained by scan- 
ning the gate voltage is believed to reflect the spatial and 
energetic distribution of microscopic resistances, it is es- 
pecially interesting to see how it is influenced by the field 
effect anomaly and its slow relaxation. We published pre- 
liminary results in Ref. [6, 38J and similar investigations 
were reported very recently on L^Oa-^ films [39 . 

In Sect. [IT] we will present samples elaboration and 



measurement techniques. Sect. Ill is devoted to the study 



of reproducible conductance fluctuations in samples of 
micrometer size. The data are shown to be in agreement 
with a percolation model applied to strongly inhomoge- 
neous media. In Sect. |IV| the conductance fluctuations 
are studied in relation to the anomalous field effect (the 
conductance dip) and slow relaxation phenomena. The 
conductance dip formation and its slow response to gate 
voltage changes are shown to have no influence on the 
fluctuations pattern. Conversely, the existence of con- 
ductance fluctuations does not modify the conductance 
dip parameters. 



II. SAMPLES AND EXPERIMENTALS 

A. Samples elaboration 

The samples used in this experiment are MOSFET de- 
vices in which the conducting channels are made of insu- 
lating granular Al thin films. To obtain them, granular 
Al is evaporated on top of a heavily doped Si wafer (the 
"gate") covered with a 100 nm thick Si02 insulating layer 
(the "gate insulator"). The evaporation was done either 
through a hand made mask (HM samples) or through a 
resist mask patterned by classical electronic lithography 
techniques (EL samples, see Fig. [TJ. In this last case, 
20 nm thick Au contact pads are evaporated prior to re- 
sist mask fabrication, whereas for hand made masks, Al 
contact pads could be evaporated on top of granular Al 
films without opening the vacuum chamber of the evap- 
orator. 

Granular Al films are deposited by e-beam evapora- 
tion of Al under a controlled pressure of oxygen. We fo- 
cus here on samples lying on the insulating side but still 
close to the metal-insulator transition. The electrical re- 
sistance per square Ro at 4 K is tuned between a few MO 
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FIG. 1: Scanning electron microscope picture of a granular 
Al channel made by electron lithography. The line is 180 nm 
wide and 820 nm long. 



to a few Gil by slightly increasing the oxygen pressure. 
Typical parameters are an Al evaporation rate of 2A/s 
and an oxygen pressure around 2.10~ 5 mbar. The film 
thickness is 20 nm for hand made samples and 40 nm for 
samples made by electronic lithography. The microstruc- 
ture of the films was already discussed in Ref [5] . They 
are believed to consist of Al grains with typical diameters 
of 5 nm separated by thin insulating A10 K barriers. 



ments. The measured conductances were very close (the 
~ 1% difference is of the order of the instruments pre- 
cision) and the gate voltage fluctuations patterns to be 
discussed in this paper were practically indistinguishable 
(see Fig. [2). 
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FIG. 2: Normalized conductance measured as a function of 
the gate voltage V 9 for a channel 20/im long and 30/^m wide 
(20/xm x 30/j.m) in a 2-contacts (triangles, average of 20 V 9 
sweeps) and 4-contacts (circles, average of 5 V 9 sweeps) con- 
figurations. HM sample, Rn = 2 Mf2 , T = 4.2 K. 



B. Electrical measurements 



The electrical resistance of the granular Al channels 
was measured in a two contacts configuration. The 
sample was DC or AC voltage biased and the resulting 
current was measured through a home made or Femto 
DLPCA-200 current amplifier. A DC voltage was ap- 
plied to the gate and could be swept between -30 V and 
30 V. No current leak through the gate insulator was de- 
tectable in this voltage range. DC voltage sources for the 
gate and the bias were Yokogawa 7651 or AO IP SN830. 
AC measurements (frequency between 10 Hz to 200 Hz, 
depending on the sample resistance) were done with Sig- 
nal Recovery 7265 or Stanford SR810 Lock In amplifiers. 

All the measurements presented below have been done 
at low bias voltages, in the linear regime of the I-V curves: 
typical bias voltage is 50 mV for a 50 fim long channel. 

Some samples made by electronic lithography dis- 
played pronounced non linear I-V effects for unexpect- 
edly low voltages compared with hand made samples of 
similar sizes. They were not used in this study. For 
hand made samples, the absence of contacts effects in 
the measured conductance was carefully checked. First, 
I-V curves were measured for samples of different lengths 
and it was checked that they depend only on R n and the 
electric field E across the granular Al channels. Second, 
for one not too resistive sample (R a = 2Mf2) 20/ito long 
and 30/ito wide, we compared 2 and 4 contacts measure- 



Ill. CONDUCTANCE FLUCTUATIONS IN 
SMALL SAMPLES 

A. A reproducible conductance fluctuations 
pattern 

When the samples size falls below sa 100/im, clear and 
well resolved fluctuations of the electrical conductance 
are visible as a function of the gate voltage at 4 K. The 
typical relative amplitude is 0.1% for a 50/ito x 50/im 
sample. 

The fluctuations pattern does not depend on the mea- 
surement conditions. It is independent of the gate voltage 
sweep parameters (sweep rate between 8 s/V to 160 s/V, 
sweep direction, etc.), of the instruments used (gate and 
bias voltage sources, current amplifiers) and of the work- 
ing frequency (DC or AC measurements). Indeed, in ad- 
dition to the test of Fig [2] the dependence of the fluc- 
tuations amplitude with the square root of the granular 
Al channel area (see Fig. |8| makes up an indirect confir- 
mation that the pattern actually comes from the channel 
itself and not from the contacts. 

As we can see in Fig. [3] the fluctuations pattern re- 
mains stable for one day after a quench at 4.2 K. Some 
noise is always superimposed on the fluctuations pat- 
tern that can be reduced by averaging successive sweeps. 
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FIG. 3: Normalized conductance measured over one day for 
a channel with a constriction of 340nm x 230nm (EL sample, 
Ra = 3 MQ , T — 4.2 K). The first gate voltage sweep (curve 
labeled t = 0) was started less than one minute after a quench 
from 10 K down to 4.2 K. The other curves have been shifted 
for clarity. 



A careful check of the long-term stability often reveals 
a small decrease of the correlation coefficient C be- 
tween successively measured patterns @0], reflecting mi- 
nor changes in the details of the fluctuations. If the sam- 
ple is warmed up to room temperature and quenched 
again, the pattern is completely different (see Fig. |4^). 
A more detailed study indicates that the temperature has 
to be raised above ~ 30 K in order to change significantly 
the fluctuations (see Fig. |I]d). 

The distribution of the conductance values G is 
roughly described by a Gaussian function. In disordered 
insulators and when the conductance fluctuations are 
large, asymmetries in the distribution of log G could give 
important information about the geometry of the criti- 
cal resistance network [411 142 j . For example, a ID chain 
will result in a tail towards low conductance values. But 
even in the sample with a narrow line 100 nm wide, we 
couldn't find any systematic asymmetry in the conduc- 
tance distribution. This is probably due to the fact that 
our fluctuations are small (1% or less of the conductance 
value ) . 

We characterized the fluctuations pattern by its root 
mean square (rms) amplitude tr (standard deviation of 
{G(V g )i} data). The noise contribution could be sub- 
tracted either by working out the average of many scans 
or by estimating its rms contribution (the variance of 
the difference between two scans gives an estimate of 
twice the variance of the noise). In order for a to be 
a well defined number, we used V 9 scans with at least 
one point every 0.5 V and over a range of 40 — 60 V. 
For a given channel, a is found to vary by about 10% 
between different coolings. Fourier Transforms of the 
fluctuations pattern does not reveal any V g periodicity 
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FIG. 4: (a) Normalized conductance for a channel with a con- 
striction of 340nm x 230nm (EL sample, R D = 3 Mfi). The 
upper curve corresponds to the fluctuations pattern at 4.2 K 
after a first cool down. The lower curve (shifted for clarity) 
corresponds to the new pattern obtained after the sample was 
warmed up to room temperature and cooled again at 4.2 K. 
The correlation coefficient C between the two curves is 0.19 
while it is around 0.92 between successive sweeps at 4.2 K. 
(b) Normalized conductance for a channel with a constriction 
of 1.8/xm x 130nm (EL sample, Ra = 3 Mfi). From top to 
bottom: the initial fluctuations pattern at 4.2 K; the 4.2 K 
fluctuations patterns after 10 mn excursions at 10 K, 20 K, 
30 K and 40 K. The correlation coefficients C indicated on 
the figure have been calculated between 4.2 K sweeps taken 
before and after the temperature excursions. The curves have 
been shifted for clarity. 



and the autocorrelation function falls down to half its 
value for V g = 2 — 5 V. A zoom of the fluctuations on 
a reduced V g scale is shown in Fig. [5j We can see that 
well-reproducible structures of only 0.5 V wide are also 
present. 

The pattern is also independent of the equilibrium 
gate voltage which is maintained between the sweeps (see 
Fig. |6|. We will come back in more details to this obser- 



5 



3,5 1 0" 9 i _ 1 | 1 1 1 | 1 1 1 | 1 1 1 i 1 1 1 | 1 




3.44 10" 9 r 1 1 1 1 1 1 11 



-2 2 4 6 



FIG. 5: Conductance for a channel of x 50/im (HM 

sample, Ra = 610 MQ, T = 4.2 K). The upper curve (trian- 
gles) is the average of 12 successive sweeps, and the bottom 
one (circles), shifted for clarity, the average of 9 sweeps. Small 
but well-reproducible structures are visible on a gate voltage 
range of 0.5 V. 
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FIG. 6: Normalized conductance for a channel with a constric- 
tion of 1.8/Ltm x 130nm (EL sample, R a = 3 Mf2, T = 4.2 K). 
The upper curve is the average of 9 sweeps taken after a long- 
time equilibrium under V g — V. The middle curve is the av- 
erage of 13 sweeps measured after the equilibrium gate voltage 
was changed to 10 V. And the bottom curve is the average of 
7 sweeps taken after the equilibrium gate voltage was changed 
to -10 V. The duration of each sweep was 240 s, with 1 hour 
between two successive sweeps. Middle and bottom curves 
have been shifted for clarity. 



B. General picture of the electrical conduction in 
our granular Al films 

We now have to wonder how to understand these con- 
ductance fluctuations and the features mentioned above. 
In particular, how could they be compared to the conduc- 
tance fluctuations observed in "homogeneous" disordered 
insulators, like indium oxide? 

As it was mentioned in Sec. |TTJ our insulating granular 
Al thin films consist of Al grains about 5nm in diame- 
ter separated by AlO^ tunneling barriers. One impor- 
tant energy scale is the charging energy of an individual 
grain. Assuming a tunneling barrier thickness of 1 nm 
and a dielectric constant close to that of AI2O3, we get 
a charging energy of about 150 K. It is much larger than 
the thermal energy at 4 K and Coulomb blockade effects 
are thus important at low temperature in our system. 
The room temperature Rn of the films spans from 30 
to 80 k£! (HM samples, 20nm thick). It should give a 
lower bound value for the dominant microscopic tunnel 
resistances [43 . Even if such values are not much larger 
than the quantum resistance h/e 2 ~ 26 kO, we will ne- 
glect quantum fluctuations of charge on the grains in the 
following discussion. 

In a perfectly regular 2D-array of metallic grains (i.e. 
an array with unique values of inter-grains capacitances 
C, gate-to-grains capacitances C 9 and tunnel resistances 
between adjacent grains), a global gate voltage V 9 in- 
duces the same polarization charge in each grain Q g = 



CgVg. A naive picture suggests that at low temperature, 
such a V g change gives rise to a large and periodic mod- 
ulation of the array conductance with a period of e/C g , 
similar to what is observed in Single Electron Transis- 
tors or SETs. This is clearly not what we are observing 
here: relative conductance fluctuations are of only few 
percent or smaller. But as we recall now, this naive pic- 
ture neglects important sources of potential disorder in 
the system [TIE]. 

The first source of disorder is the inevitable presence 
of grains size and tunnel barriers thickness distributions, 
which unfortunately we cannot measure easily. They 
would result in distributions of charging energies, gate-to- 
grains capacitances C g and tunnel resistances. But even 
if a wide distribution of C g will smooth the Coulomb 
oscillations for gate voltage larger than w e/C g , a 
large gate voltage conductance modulation should still 
be present at V g — (grains in phase) which is again not 
observed in our films. 

Another important source of random potential comes 
from the presence of charged impurities or defects 
trapped in the dielectric environment of the grains (tun- 
neling barrier, substrate and natural oxide layer cover- 
ing the film). Such trapped charges polarize the metal- 
lic grains and induce so called offset charges. The off- 
set charges Qo are not quantized and could be larger 
than e. They can be static or dynamic. When a metal- 
lic island is connected to bulk electrodes with two tun- 
nel junctions (SET), the current through the device at 
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low bias (in the Coulomb blockade regime) was shown 
to be very sensitive to tiny changes in the island po- 
tential corresponding to offset charges of 10 _3 e or be- 
low [HI SB]. The existence and the dynamics of off- 
set charges were clearly identified in metal-based transis- 
tors made by electronic lithography techniques (typical 
island size is 1 /zm) @5J S3 EH HO M EE], in scan- 
ning tunneling spectroscopy measurements on individual 
and oxidized metallic grains of nanometer size [S3 [53] 
and in planar metal-based tunnel junctions [54 1 55] . In 
our experiments, the granular Al films are deposited on 
a thermally grown SiC>2 layer 100 nm thick. We think 
that, in agreement with results on Al oxidation |56j . our 
AlOa; dielectric layer around the metallic grains is not 
stoichiometric with AI2O3 and contains a large concen- 
tration of oxygen vacancies. This layer certainly consti- 
tutes the main source of potential disorder in our system. 
The role of interface states in the electronic properties of 
Al-AlOj; Al tunnel junctions was also emphasized in Ref. 
[21] . Capacitive studies on standard Si-SiC>2 wafers cov- 
ered with an electron-gun evaporated Al layer, indicate 
typical charge densities about 10 12 e.cm -2 , i.e. about one 
charge every 10 nm (mainly interface states) |57j . 

One last source of disorder for the potential of the 
grains was considered in Ref. j58j [59] . When the grains 
are very small, surface effects (irregular shapes and sizes) 
contribute to random changes in the grains potential that 
are much larger than the mean energy level spacing A, 
and that could be even larger than the charging energy 
Eq. For metallic grains of 5 nm in diameter, Cuevas et 
al. [55J estimate that this effect could by itself ionize 
more than half of the grains. In our case, by taking the 
approximate formula A = l/n(Ep)d 3 (d is the diameter 
of the grains) and the density of states of bulk Al, we get 
that A is about 4 K. 

All these effects together result in a large distribu- 
tion of offset charges that could exceed many times the 
electron charge e. Such large offset charges distribution 
is usually unstable with respect to single electron tun- 
neling between the grains. They are (partially) com- 
pensated by quantized electrons hops from other grains 
and/or from the electrodes. Theoretical studies on regu- 
lar arrays of tunnel junctions have shown that the offset 
charges configuration that minimizes the electrostatic en- 
ergy depends on the ratio C/C g [BUIGS]- In ID and when 
C g <C C, a charge placed on one island is screened on a 
distance (C/Cg) 1 / 2 called the soliton length (in units of 
the array step-size) [62]. When C <C C g (short screen- 
ing length limit), it is usually assumed that the system 
ends with a random and uniform distribution of offset 
charges between —e/2 and e/2 [THJES]- But in the oppo- 
site limit (C 3> C g ), numerical simulations on ID arrays 
have found that starting from a random offset charge 
distribution, the interactions between electrons on dis- 
tant grains induce energy favorable electron tunneling 
between the grains that smooth the potential [60]. Simi- 
lar simulations were also done on 2D-arrays smaller than 
the soliton length [51]. After electrostatic energy mini- 



mization, the density of states of single electron addition 
energies has a Coulomb gap equivalent to what is ob- 
served in homogeneous disordered systems (but with a 
different energy dependence). Moreover, the linear con- 
ductance at very low bias voltages (eV <C ksT) is ther- 
mally activated. The activation energy is about 0.1e 2 /C 
and changes in the offset charges configuration induce 
activation energy fluctuations about l/3rd of this value 
(these quantities are almost independent of the size of 
the system). In our granular Al thin films, a rough esti- 
mate gives C = 2.10~ 18 F and C g = 7.10~ 21 F (soliton 
length (C/Cg) 1 / 2 — 20). Thus, electronic correlations 
between the grains may play an important role, even if 
the soliton length (about 100 nm) is still smaller than the 
system size. It is worth noting that the understanding 
of the electronic properties of granular systems is still an 
active theoretical field [53] in spite of decades of research. 

For simplicity, we will consider now that we are let 
with a random potential disorder with a homogeneous 
distribution of offset charges on the different grains (we 
neglect electronic correlations). Such a distribution nat- 
urally explains why large Coulomb oscillations are not 
observed as a function of V g in macroscopic samples [44] . 
It also implies (in the regime kgT <C Ec) an exponen- 
tially large distribution of hopping probabilities between 
the grains (a distribution in the barrier thicknesses will 
have the same effect) [65], Therefore, percolation theory 
results obtained for strongly inhomogencous media [26 
should also describe the conductance of our granular Al 
films. According to this theory, our V g fluctuations come 
from an incomplete self-averaging of the conductance in 
small size samples (" mesoscopic" fluctuations). Like in 
homogeneous hopping systems, the conductance is dom- 
inated by a small number of critical resistances (which 
are here single electron transistors, possibly asymmetric) 
with a typical distance between them given by the cor- 
relation length of the critical network (see below). Our 
stable fluctuations pattern reflects a specific distribution 
of offset charges on the grains, which we call below the 
background charge distribution. A gate voltage sweep 
shifts all the offset charges and changes the conductance 
and/or the position of the dominant single electron tran- 
sistors. 

Annealing the system to room temperature resets the 
background charge distribution and completely modify 
the fluctuations pattern as seen in Fig. |4^ (a direct ev- 
idence of annealing effects in individual grains was seen 
in Ref. |66j). Our results for annealing at intermediate 
temperatures (Fig. [4jj) may indicate that the background 
charge distribution is thermally activated. We often ob- 
serve at 4K a small decrease of the correlation coefficient 
between V g scans as time has elapsed, which may corre- 
spond to long-term drifts of offset charges at low T. 

The absence of periodicity in G(V S ) fluctuations could 
be simply explained by a distribution in C g values. Sig- 
nificant conductance modulations occur on a scale as 
small as V 9 = 2 V, 10 times smaller than the typical 
scale e/ <C g >. It corresponds to an average change in the 
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charge per grain of only O.le. If such a strong sensitivity 
to V g changes could also come from a large C g distribu- 
tion, we believe that the interactions between electrons 
on distant grains may play an important role. Indeed, 
numerical simulations on ID regular arrays have shown 
that if the soliton length is large, the threshold voltage is 
sensitive to gate voltage changes that represent a fraction 
of electrons per grains [60 . Moreover, the modulation 
pattern of the conductance is strongly dependent on the 
actual background charge distribution. Unfortunately, 
measurements of 2D artificial junction arrays made by 
electronic lithography are not directly comparable to our 
work [H7J EH]- They do not study V 9 conductance fluctu- 
ations in the linear regime but above the threshold volt- 
age (non ohmic regime) and the interactions are usually 
short-range (C/C g below or close to 1). 

C. Size and R n dependence of the conductance 
fluctuations 

In order to measure how the fluctuations amplitude 
depends on the samples size and resistance, we have pro- 
duced a large set of HM samples. Sizes are distributed 
between 15 /mi and 1 mm, and Rn at 4.2 K between 1 Mfi 
and 1 Gil. Unless explicitly mentioned, all the following 
measurements have been done at 4.2 K. We will show 
that, in agreement with the previous discussion, most 
of our results can be interpreted within a general perco- 
lation model applied to a system with an exponentially 
large distribution of microscopic resistances. 

1. Size dependence of the conductance fluctuations 

For samples having similar R n , the fluctuations ampli- 
tude is strongly dependent on the sample size, the larger 
the channel area, the smaller the amplitude (see Fig. [7]). 
Indeed, reproducible fluctuations are still measurable on 
samples of mm size but the rms relative amplitude falls 
below 0.01%. 

Relative rms amplitudes are gathered on Fig. [8] for 
three sets of samples with respectively Ra about 5 Mf2 , 
50 Mfi and 500 Mil . The relative rms amplitude is found 
to decrease as the invert of the square root of the channel 
area. This dependence indicates that the observed fluc- 
tuations result from the sum of statistically independent 
microscopic fluctuations. The length scale L at which 
(G- < G >)/ < G >= SG/ < G >= 1 is given by 

i fIo ««(2D) (^) = (^)" 2 (1) 

i f X » i( 3D, (^k) = ( L 4f < 2 > 

where t and S are respectively the thickness and the 
area of the film. Assuming a 2D regime, we get from the 
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FIG. 7: Normalized conductance for channels of different size 
(HM, R a = 200 MQ, - 1 GO). Channels with fluctuations 
amplitudes in decreasing order: 20/im x 15^im, 50fim x 35/tm, 
llO^im x 105/im, 1.1mm x 50/tm, 1.1mm x 870/im.. 

data of Fig. [8] L = 30 nm for R n = 5 Mf2 and 60 nm 
for R a = 500 MO. An interpretation of the length Lo 
comes from the percolation theory applied to strongly 
inhomogeneous media Let's suppose that we have 

an exponentially wide range of resistances Rij between 
the sites i and j of an array, i.e. Rij — exp(fy) with an 
uniform distribution of £jj between < £y < A£. Then, 
if A£ ^> 1, the resistance R of the array is that of the 
critical resistance subnetwork and is given by 

i? = i? exp(£ c ) (3) 

The exponential factor £c is the smallest ^ that first 
gives percolation when only sites such as ^ < £c are 
connected. The critical resistance subnetwork includes 
all the sites such as £ < £c + 1 (sites with larger 
have a negligible contribution) . It is characterized by its 
correlation radius or homogeneity length Lo- In other 
words, one unavoidable microscopic resistance Rjj of the 
order of the macroscopic one is encountered on a typical 
length scale Lo. Lo is given by 

Lo « W (4) 

v = 0.9 at 3D, v = 4/3 at 2D and I is the microscopic 
length of the problem (array step size if the hopping is be- 
tween nearest neighbors or hopping length in the variable 
range hopping regime). Changing V g shifts the Fermi 
level of the system and therefore the values of the mi- 
croscopic resistances in the system. This gives rise to 
conductance fluctuations SG/ < G >= 1 at the scale L 
and to modifications of the critical resistance network. 
We will see below whether the condition A£ ^> 1 neces- 
sary for the validity of Eq. [3] and Eq. [4] are fulfilled in our 
samples. 
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FIG. 8: Relative rms amplitudes of the conductance fluctu- 
ations as a function of the channels areas S (HM samples, 
T = 4.2 K)). Three different set of samples with similar R n 
are represented: Rn ~ 5 Mf2 (circles), Ra ~ 50 Mil (squares) 
and Ra ~ 500 Mil (triangles). The straight lines correspond 
to the square root dependence of Eq. [7J and [5] (see text for 
details). 



The Lo values extracted from data of Fig. [8] are larger 
than the film thickness (20 nm for HM samples), con- 
firming a 2D- like regime. Our Lo values are signifi- 
cantly smaller than the homogeneity lengths obtained in 
Ir^Oa-z thin films (Lo ~ 300 nm) (39]. A study of the 
fluctuations amplitude dependence on the film thickness 
(especially around Lo) would be very interesting in order 
to confirm our estimates. The role of the electronic inter- 
actions between distant grains on the percolation picture 
developed above, especially for length scales smaller than 
the soliton length, should also be clarified. 
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FIG. 9: Relative rms amplitudes of the conductance fluctua- 
tions as a function of Rn .The amplitudes correspond to chan- 
nels area S of 50 \im x 50 jj,m: they have been deduced from 
data on channels with different sizes and assuming the square 
root dependence of Eq. [7J The dispersion of the data gives an 
idea of the rms amplitude uncertainty. 
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2. Ra dependence of the conductance fluctuations 



1/T(K- n ) 



Fig. [8] also emphasizes an increase of the conductance 
fluctuations amplitude as a function of R n for a given 
channel area. The increase is weak: the relative rms 
amplitude is multiplied by only 2 when R n rises from 
5 MO to 500 Mtt (see also Fig. [9). 

Below ~ 10 K, the resistance temperature dependence 



is close to an activated law (see Fig. 10 1 



R = R exp(E /k B T) 



(5) 



According to the percolation model described previ- 
ously, 



£c = E Q /k B T » AC 
Eq. |4] and [6] imply that 

L « l{E /k B T) v 



(6) 



(7) 



FIG. 10: R D as a function of 1/T between 4.2 K and 300 K 
for three HM samples. The low T part (4.2 K-10 K) could be 
well described by an activation law. 



Since for a given channel area the fluctuations amplitude 
is proportional to Lq (2D-regime) we expect at fixed T: 



5G 

<G> 



(8) 



In Fig.[TT] the rms fluctuations amplitudes are plotted as 
a function of the activation energies Eo extracted from 
the fits between 4.2 K and 10 K. The data are in good 
agreement with Eq. [8] taking v = 4/3. The activation 
energies range between 16 K and 35 K for HM samples 
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with R a = 3 Mil and 1 Gil respectively. This justify the 
use of Eq. [3] and [I] since A£ = £c/%c — EqKUbTxc) is 
much larger than 1: A£ > 4 for R n = 3 Mil and A£ > 8 
for Rn — 1 GO (xc, the critical probability, depends on 
the array type but is always smaller than 1). 



decreases as T is increased (see Fig. |L2ft). The decrease 
is well described by the power law Eq. [9] but the expo- 
nents a are above the expected ones (cv=2 and 1.7 for 
respectively the HM and the EL sample, see Fig. 12 d) 



0.14 




1.Q2 



FIG. 11: Relative rms amplitudes of the conductance fluctu- 
ations (area S of 50/xm x 50/im) as a function of E^ 3 . The 
straight line corresponds to Eq. [8] 

Up to now, the resistance T dependence of our insu- 
lating granular Al thin films is still not fully understood 
in the whole T range and interactions between electrons 
may play a role. The activated law of Eq. [5] is a reason- 
able approximation of the low T behavior. Strictly speak- 
ing the standard percolation model developed above ap- 
plies to non-interacting electrons. It is anyway outstand- 
ing that the above analysis gives a coherent picture of 
the data. 
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D. Temperature dependence of the conductance 
fluctuations 

We have also measured the conductance fluctuations 
as a function of T above 4K for two samples: one HM 
sample with R a = 1 Gil (Eq ~ 35 K) and one EL sam- 
ple with R n = 3 Mil (E Q ~ 20 K). A higher T implies a 
resistance exponentially smaller. In terms of the percola- 
tion theory scenario, the correlation length is smaller, the 
number of dominant microscopic resistances increases in 
the system which finally reduces the relative fluctuations 
amplitude. According to Eq. [7] we expect: 



SG 

<G> 



cc (l/TY 



(9) 



With a = v = 4/3 at 2D, a = (3/2) 1/ = 1.35 at 3D. 

The results for the two samples are presented in Fig. [12] 
and Fig. |13| The agreement with the percolation theory 
predictions is only qualitative. The relative amplitude 



FIG. 12: (a) Relative conductance as a function of V 9 for 
different temperatures. The sample is a channel constriction 
of 1.8fj,m x 130nm (EL, Rn = 3 MQ,). By order of decreasing 
fluctuations amplitude: T = 5.2 K, 10 K, 20 K, 30 K and 40 K. 
(b) Relative rms fluctuations amplitude as a function of T. 
Squares: HM sample with channel size 115/imx 105/xm (Rn = 
1 Gf2); straight line: power law Eq. [9] with a = 2.0. C ircle s: 
EL sample with Rn = 3 MQ (same channel as in Fig. |12^ ) ; 
straight line: power law Eq. [9] with a = 1.7. The full and 
empty circles represent two set of measurements performed 
on two different cryogenic sticks. 



As illustrated on Fig. 13 1 and b, the absolute rms 
fluctuations amplitude first increases and then decreases 
slightly above 10-20 K. In Fig. [ilfc, it is seen for the EL 
sample that the positions of the main structures are un- 
changed below 20 K (the autocorrelation function width 
at half value is indeed constant in that T range) and 
that they are smoothed significantly above. If we still 
don't know how to explain the observed T dependencies 
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of Fig. 12 and Fig. 13 the extension of the measurements 
down to lower T may bring interesting information on 
this problem. 
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FIG. 13: (a) Absolute conductance variations G— <G> as a 
function of V 9 for different temperatures. Same sample as 



in Fig. [12k, From top to bottom: T = 5.2 K, 10 K, 20 K, 
30 K, 40 K. The lower curves have been shifted for clarity, 
(b) Absolute rms fluctuations amplitude as a function of T. 
Same samples as in Fig. |12b . 



IV. THE ANOMALOUS FIELD EFFECT IN 
PRESENCE OF THE CONDUCTANCE 
FLUCTUATIONS 



One main goal of these experiments was to study the 
anomalous field effect and its slow relaxations in pres- 
ence of the conductance fluctuations. Slow conductance 
relaxation phenomena in our macroscopic (with negligi- 
ble fluctuations) granular Al films have been described in 
details in Ref. [5J [5] . 



A. Time evolution of the conductance dip 

In macroscopic samples maintained under a fixed gate 
voltage Vpetf at 4.2 K, a fast gate voltage sweep reveals a 
conductance dip in the G(V ff ) curves, symmetric and cen- 
tered on V ffeg . The dip amplitude was found to increase 
as a logarithm of the time elapsed since V g = V geq . The 
formation of a conductance dip after a quench at 4.2K 
under V geq = V is illustrated in Fig. [lih. The same 



protocol has been followed in Fig. [T4p but on a smaller 
sample (50/wn x 40 [im channel). The conductance dip is 
still present but superimposed on the conductance fluctu- 
ations pattern described previously. The two effects are 
here of equal magnitudes. The striking feature is that the 
fluctuations pattern is not affected by the digging of the 
conductance dip, neither for V g values within the con- 
ductance dip range (±3 V around V seg ) nor out of this 
range. Actually, long-term changes are also visible in the 
fluctuations of Fig. [l4p , but they are weak and similar to 
what is sometimes seen after a quench in samples where 
the conductance dip is not visible. Since the fluctuations 
and the noise amplitude increases when the sample area 
is reduced (which is not the case for the dip amplitude, 
see Sect. IV B| , the conductance dip is not visible in too 
small samples (channels size below ~ 10[im). This was 
for example the case for the EL sample pattern of Fig. [6] 
We have also tested a second procedure. Once the 
sample has been kept for a long time (days) under V geq \, 
the gate voltage is changed to a different value V geq 2- 
The results are illustrated for a macroscopic sample in 



Fig. 15 1: a new dip is formed at V ffe <j2 while the old dip 
* v gegi is slowly erased. In microscopic samples (see 



Fig. 15 d), the conductance dip behaves similarly and the 
fluctuations pattern is essentially unchanged under the 
displacement of the conductance dip. 

This is even clearer in Fig. |16h , and b, where similar 
data are presented for a microscopic sample but on a re- 
duced V g range around V geq 2- By subtracting the curve 
at t = (measured under V geq i) from the other curves 
(measured under V geq 2)j the stable fluctuations pattern 
is eliminated and we are let with a smooth and symmet- 
rical conductance dip (see Fig. 16 a). Thus, like for the 
previous procedure, the fluctuations pattern is not sig- 
nificantly modified by the formation or the erasing of a 
conductance dip, for gate voltages inside or outside the 
conductance dip range. We have carefully checked this 
point in many samples having different Rn values and 
down to an area of 20/im x Ibfim. 



B. Parameters of the conductance dip: 
amplitude 



width and 



We have compiled conductance dip amplitude data 
measured on samples of different sizes (lowest area 
20/im x 15/im) and different Rn values. Since the dip 
is an off-equilibrium property, it is crucial to use exactly 
the same procedure for all the measurements. In Fig. |17| 
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FIG. 14: (a) Normalized conductance for a "macroscopic" 
sample (HM, Ra = 30 MQ, 100/im x 2mm). The sample was 
quenched from room temperature to 4.2 K at t — under a 
fixed gate voltage of V geq = V. Scans of 240 s long were taken 
every 30 mn. A symmetric dip centered on V geq increases in 
amplitude with time. The lower curves have been shifted for 
clarity, (b) Normalized conductance for a "microscopic" sam- 
ple (HM, Ra = 3.3 MQ, channel size 50pm x 40^m). Like 
in (a), the sample was quenched from room temperature to 
4.2 K at t — under a fixed gate voltage of V geq — V. Scans 
of 240 s long were taken every 2 hours. Reproducible conduc- 
tance fluctuations are superimposed on the conductance dip. 
The lower curves have been shifted for clarity. 



FIG. 15: (a) Normalized conductance for a "macroscopic" 
sample (same sample as in Fig. 14 1). The sample was main- 



tained 24h under V„ 



-7.5 



V and the gate voltage was 



changed to +7.5 V after the scan labeled t = s. 80 mn 
later, a new dip is present at 7.5 V while the amplitude of 
the dip at -7.5 V has been reduced. The lower curve has been 
shifted for clarity, (b) Normalized con duc tance for a "micro- 
scopic" sample (same sample as in Fig. 14 a). The sample was 
maintained 4 days under V geq i = V and the gate voltage 
was changed to V geq 2 = +10 V after the scan labeled t = 0. 
A new dip which amplitude increases with time is formed at 
+ 10 V while the "old" dip at V is erased with time. Scans 
are 80 s long, with 400 s between two scans. The lower curves 
have been shifted for clarity. 



the parameters are the following: the samples are first 
maintained under a fixed V ge(? i for at least 24 h. The 
gate voltage is then changed to V ge q2 for 1 hour, creat- 
ing a new dip, and a gate voltage scan 20 V wide starting 
from V geq 2 is recorded (81 points, 4 s/point). The differ- 
ence between the conductance out of the dip range and 
the conductance at V ge q2 gives the new dip amplitude. 
For microscopic samples (where significant fluctuations 
are present), a reference scan taken before the change 
to V ge9 2 is subtracted. The results are independent of 



the Vgeqi and V ge q2 values. Fig. 17 clearly demonstrates 



that, in spite of a larger data dispersion in microscopic 
samples (due essentially to a larger noise level), the con- 
ductance dip amplitude depends on R n and not on the 
sample area in contrary to the fluctuations amplitude. 
Moreover, the R n dependence of the dip amplitude is 
faster than that of the conductance fluctuations. Indeed, 
its relative value is multiplied by 10 between 2Mf2 and 
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FIG. 16: (a) Normalized conductance for a "microscopic" 



sample (same sample as in Fig. 15 d). The sample was main- 
tained three days under V geq i = V before the gate voltage 
was changed to V geq 2 = 10 V after the scan labeled t = 0. 
The V 9 values of the sweeps focus on the dip range, (b) Dif- 
ference between the G(V g ) curves at time t > (measured 
under V geq 2) and at t = (measured under V geq i). 



IGil , instead of 3 for the fluctuations amplitude (see 
Fig. [9). 

The conductance dip width is also the same in the 
different samples we have measured, independently of the 
conductance fluctuations amplitudes and the R n values 
(see Fig. 18 1. As it was already mentioned in Ref. [BJ, 
the width of the dip in our granular Al films is fixed by 
the temperature of the measurement and is independent 



of R n values and V„ sweep parameters. 



C. Time relaxation of the conductance 

Finally, we have tested if the conductance relaxation 
laws were modified by the presence of conductance fluc- 
tuations. One sample is maintained under a fixed gate 




R (O) 
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FIG. 17: Relative amplitude of the conductance dip for dif- 
ferent HM samples. The amplitudes correspond to that 
of dips dug for one hour at 4.2 K (see text for details). 
The microscopic (circles) and macroscopic (squares) samples 
have respectively areas in the range 1CF 10 — 1CF 8 m 2 and 
10-8 -10~ 6 m 2 . 



voltage during a few days and we follow the conductance 
response to a sudden change of V g . In macroscopic sam- 
ples [5J, after a fast increase, the conductance decreases 
as a logarithm of the time elapsed since V ff was changed. 
In presence of conductance fluctuations, the decrease is 
noisier but is still described by a logarithm of time (see 
Fig. pi. 



We have also tested the so called two dips experiment 
[3]. In this protocol, the sample is maintained for a long 
time under V geq i, a new dip is formed at V ge(? 2 during 
a time t w and the evolution of G(V g ) curves is followed 
when the gate voltage is switched back to Vgeqi- At least 
when | Vg eq i — Vgeq2 | is not too large and in macroscopic 
samples, the erasing of the V geq 2 dip amplitude was found 
to be a universal function of t/t w (simple aging) 3,6 . 
Once again, we have observed the same relaxation curves 
in macroscopic and microscopic samples (see Fig. [20| , 
whatever the conductance fluctuations amplitude. 



D. Discussion 

Our measurements demonstrate the following impor- 
tant features: 

- the glassy conductance dip is not significantly al- 
tered by the samples size down to 20 /im. This 
suggests that the conductance dip doesn't result 
from an incomplete self-averaging of random micro- 
scopic conductance modulations. Its shape and am- 
plitude are indeed independent of the sample size 
even when the fluctuations amplitude is of the same 
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FIG. 18: Shapes of the conductance dip for HM samples of 
different sizes. The stable fluctuations have been removed 
by subtracting a reference G re /(V 9 ) scan taken before the 
formation of the dip. The fluctuations that remain (larger 
in smaller samples) are due to the short time noise level of 
the scans. Channels sizes are: 1.1mm x 50[im (R = 60 MS7 ), 
50/imx40/^m (R = 6.6 Mtt ), and 20m^x 15/im (R = 240 MQ 
). The conductance was arbitrarily set to in its minimum 
and to 1 out of the dip range. 



FIG. 20: Two dips experiment measurement in a macro- 
scopic and a microscopic samples. We plot how the V geq 2 
dip is erased with time once V g is switched back to V geq i 
[AG =< G > — G(V geq 2), see text for details]. Macroscopic 
sample: HM, R a = 30 MQ , channel size 100/im x 2mm, 
t w = 3600 s. Microscopic sample: HM, Ra ~ 3.3 Mf2, channel 
size 50/xm x 40/im, t w = 10000 s. The microscopic relaxation 
curve has been shifted and normalized on the macroscopic one 
for comparison. 
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FIG. 19: Relative amplitude of a new conductance dip a 
time t after the gate voltage is changed from V ge9 i to V ge ,2 
[AG =< G > — G(V g eq2)] for a macroscopic (lower curve) and 
a microscopic (upper curve) samples. Macroscopic sample: 
HM, Ra = 30MS1 , channel size 100/im x 2mm. Microscopic 
sample: HM, Ra = 450M , channel size 45/im x 50/^m. The 
V S e 9 2 dip amplitudes increase as a logarithm of time (straight 
line). The values are different for the two samples according 
to their different Ra values. 



order of magnitude. The uncertainty in the con- 
ductance dip parameters increases when the sam- 
ple size is reduced, but this increase comes mainly 
from a higher noise level 

- the conductance fluctuations and the dip seem to be 
independent from each other. Indeed, the fluctua- 
tions pattern is not affected by the formation of the 
conductance dip, and conversely, the conductance 
dip is formed whatever the conductance fluctua- 
tions amplitude. If we assume that the fluctuations 
pattern reflects the specific potential landscape of 
the grains, then we conclude that the formation of 
the conductance dip does not induce large changes 
in this potential landscape 

- the conductance fluctuations and the dip have dif- 
ferent dynamics: the former are essentially static 
once the sample has been quenched at 4 K, while 
the dip slowly grows. As we have seen, a dynamics 
of the fluctuations is induced by heating the sam- 
ples, while no signature of thermal activation can 
be observed in the glassy dynamics of the dip jS] . 

The independence of the conductance dip and the con- 
ductance fluctuations was also demonstrated recently on 
Ir^Oa-z thin films down to channel sizes of 2 fj,m and 
with a maximum relative amplitude of the minimum 
about 10% [3j5]. Our results underline once again the 
strong similarities between granular Al and indium-oxide. 
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What conclusions may be drawn from these experi- 
mental observations about the origin of the glassy dip? 
In a simplified intrinsic scenario, the stable fluctuations 
pattern reflects the static distribution of charges in the di- 
electric material around the grains whereas the slow con- 
ductance decrease (dip formation at constant V g ) comes 
from slow correlated electron hops. This scenario would 
naturally explain why the two phenomena are indepen- 
dent and have different dynamics, so that our observa- 
tions are in qualitative agreement with an electron glass 
interpretation of the conductance dip. But this simplified 
picture is certainly not exact in our system. As already 
discussed in Sect. |IIIB for a regular array of tunnel junc- 
tions, trapped charges in the dielectric materials around 
the grains result in a random distribution of offset charges 
between — e/2 and e/2 only when the grains are indepen- 
dent (C <C C g ). However, this distribution does not min- 
imize the global electrostatic energy when C > C g due 
to interactions between the electrons on distant grains. 
Tunneling of single electrons between the grains permits 
a large and fast diminution of the electrostatic energy 
of the system |61| . The random potential is strongly 
modified by these fast relaxation processes: the potential 
landscape is smoothed [60J, strengthening the Coulomb 
blockade and opening a Coulomb gap in the density of 
states of single electron addition energies on the grains 
|61) . Thus the random potential (and the conductance 
fluctuations) are partly of electronic origin and cannot 
be considered decoupled from the electronic degrees of 
freedom. The conductance dip formation is attributed 
to the slower part of the electronic degrees of freedom, 
consisting of correlated electrons hops [IB]. Correlated 
electron hops are still decreasing the global electrostatic 
energy but they are not able to change significantly the 
local electrostatic order. As long as the conductance dip 
relative amplitude is small, it seems reasonable to as- 
sume that the slow correlated electron hops only induce 
small shifts in the grains potential (or possibly rare large 
shifts). For a given V g , the critical resistance subnet- 
work is therefore essentially fixed by the trapped charges 
in the dielectric material around the grains and by the 
fast single electrons hops between the metallic grains. 
In other words, we expect that the conductance dip for- 
mation affects only slightly the conductance fluctuations 
pattern. The apparent observed independence of the two 
phenomena may thus simply come from the fact that the 
conductance dip amplitude is still small. To this respect, 
similar measurements on samples having larger conduc- 
tance dip amplitudes would be of utmost interest. 

According to the extrinsic scenario, the conductance 
dip formation is the result of a slow response of the back- 
ground charges around the grains. In a continuous model, 
this is described as a slow linear polarization of the dielec- 
tric surrounding the charged grains. It could be shown 
by a subtle argument that such a polarization shifts the 
grains potentials in a way that decreases the electron 
mobility |19j . The model was initially developed for an 
electron transport perpendicular to a 2D-array of inde- 



pendent grains: starting from an uniform distribution of 
offset charges between —e/2 and +e/2, the polarization 
completely redistributes the offset charges values, the fi- 
nal effect being to reduce the proportion of grains with 
offset charges close to +e/2 or —e/2 (the "unblocked" 
grains) . Even if a detailed theory for the transport along 
a 2D- array of interacting grains is still missing, we think 
that the polarization could similarly reduce the propor- 
tion of low resistance hops in our granular films |6J . Since 
the dip amplitude is small and since it doesn't result from 
the incomplete self-averaging of a random and strongly 
fluctuating microscopic quantity, it is reasonable to as- 
sume that the part of the polarization responsible for 
the conductance dip induces either small shifts in the 
grains potentials or larger shifts but for a small num- 
ber of grains. If such shifts are triggered off by changes 
in the charge configuration of two level systems nearby, 
we believe that the second hypothesis is more realistic. 
According to Sect. |HI[ the conductance fluctuations are 
instead associated with large changes in the background 
charge distribution and consequently in the critical resis- 
tance subnetwork. We could then naturally understand 
why the rare events associated with the formation of a 
small conductance dip could not influence significantly 
the fluctuations pattern. However it is more difficult to 
explain the difference in the conductance dip and fluctu- 
ations dynamics. A more quantitative approach of the 
extrinsic scenario (e.g. via numerical simulations) would 
be very helpful to sec whether it can indeed reproduce 
the coexistence of the dip and the fluctuations. 



V. CONCLUSION 

In conclusion, we have observed at 4 K reproducible 
conductance fluctuations as a function of gate voltage in 
insulating granular Al thin films of micrometer size. Such 
fluctuations are stable over days of measurement. Their 
amplitude depends on the size and the resistance per 
square of the films in agreement with a percolation model 
applied to an exponentially wide distribution of micro- 
scopic resistances. Such wide distribution comes from 
the combined effects of a charging energy larger than the 
thermal energy and an important potential disorder of 
the grains. To the accuracy of our measurements, the 
conductance fluctuations pattern is not affected by the 
slow formation of the conductance dip. Conversely and 
down to sizes of ~ 20 fim, the conductance dip parame- 
ters are unchanged even when the two phenomena have 
the same amplitudes. This apparent independence may 
result mainly from the fact that the conductance dip am- 
plitude corresponds to small changes of the macroscopic 
conductance (no more than few percent in our case) . Our 
results may be more readily interpreted in the intrinsic 
context (electron glass origin of the conductance dip), but 
further theoretical analysis would be needed to definitely 
rule out the extrinsic scenario. 
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